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Abstract 

We analyze the neutrino Yukawa matrix by considering three constraints: the out-of-equilibrium 
condition of lepton number violating process responsible for leptogenesis, the upper bound of 
branching ratio of lepton flavor violating decay, and the prediction of large mixing angles using the 
see-saw mechanism. In a certain parametrization with bi-unitary transformation, it is shown that 
the structure which satisfies the constraints can be characterized by only seven types of Yukawa 
matrices. The constraint of the branching ratio of LFV turns out as a redundant one after applying 
other two constraints. We propose that this parametrization can be the framework in which the 
CP asymmetry of lepton number violating process can be predicted in terms of observable neutrino 
parameters at low energy , if necessary, under assumptions following from a theory with additional 
symmetries. There is an appealing model of neutrino Yukawa matrix considering the CP asymmetry 
for leptogenesis and the theoretical motivation to reduce the number of free parameters. 

PACS numbers: 13.35.Hb, 12.15.Ff, 14.60.Pq 



I. INTRODUCTION: LEPTOGENESIS 



The observed Baryon Asymmetry in Universe(BAU) can be explained with a process that 
satisfies three conditions: Baryon number(B) violation, Charge conjugation(C) violation and 
Charge conjugation and Parity(CP) violation, and for the B-violating process to be out-of- 
equilibrium |l| . Leptogenesis is a scenario for BAU in which the initial lepton asymmetry is 
recycled into a baryon asymmetry by the sphaleron process |2| . The initial lepton asymmetry 
can originate from loop-level processes involving the Yukawa couplings of Majorana neutri- 
nos [3]. The baryon asymmetry AB can be generated after sphaleron process in thermal 
equilibrium washes A(B + L) out of A(B — L) + A(B + L) , where the A(B — L) is equivalent 
to the initial lepton asymmetry — AL ^, U, 5|. The constraints on chemical potentials of 
particles in thermal equilibrium establish the relations among different asymmetries. The 
amount of baryon asymmetry is described in terms of the ratio of particle number density 
to entropy density, i.e., baryon number density with respect to comoving volume element as 

Y B = <2Z=2gL, (1) 

with entropy s, the range of which is (f| 

Y B « (0.6 - 1) x KT 10 (2) 



The baryon asymmetry is related to the lepton asymmetry |5|, |7| as 



Y B = clYb-l = — -—tYl, (3) 
a — 1 

where 

= 8N F + 4N H 
a ~ 22N f + 13Nh [ ' 

for example, a = 28/79 for the Standard Model(SM) with three generations of fermions 

and a single Higgs doublet, N F = 3,N H = 1, and a = 8/23 if N F = 3,N H = 2 as in 

super symmetric models. Thus, the amount of the original lepton asymmetry should be not 

less than ~ 1CT 9 to provide necessary amount of baryon asymmetry. 

The generation of a lepton asymmetry also requires the CP-asymmetry and out-of- 

equilibrium condition. The Y L is explicitly parameterized by two factors, e,, the size of 

CP asymmetry, and k, the dilution factor from washout process. 

Y L = {nL - n ^=K^ (5) 
s g* 



where g* ~ 110 is the number of relativistic degree of freedom. The e» is the magnitude of 



CP asymmetry in decays of heavy Majorana neutrinos 



03, 



t(u r -> i#) - r(^ -> 

where z is 1 to 3 for generation. When one generation of the right neutrinos has a mass far 
below the masses for the other generations, i.e., Mi -C M 2 ,M 3 , the ej in Eq.© reduces to 
ei from the decay of Mi, 

where 3^v is the matrix of Yukawa couplings in a weak basis while the right-handed Majorana 
mass matrix is diagonal. Such reduction to ei is possible since the asymmetry at M3 was 
washed out by the process at M? and, subsequently, the asymmetry at M 2 was again washed 
out by the process at Mi Q, Q, Q 

The k in Eq.fjSJ) is determined by solving the full Boltzmann equations. The k can be 
simply parameterized in terms of K defined as the ratio of Ti the decay width of v R \ to 
H the Hubble parameter at temperature Mi [7], where K < 1 describes processes out of 
thermal equilibrium and k < 1 describes washout effect, 

K — 7T-77, (8) 



for 10 < K = Ti/H < 10 6 , and 



— K ~ 



(9) 



2V^ 2 + 9' 

for < K ^ 10 "ij . Even if the decay is much slower than expansion rate of the 

universe, there is an unavoidable wash-out effect by k ~ 0.1 [15j, the size of which does not 
vary significantly as K varies to 10. Thus, the lepton asymmetry in Eq.fjSJ) requires a lower 
limit in CP asymmetry, to say e\ > 10~ 6 considering Yl 10~ 9 . 

It has been verified that the decay processes of heavy Majorana neutrinos can generate a 
sufficient amount of asymmetry and the asymmetry could be kep t from washout in models 

with neutrinos based on the see-saw mechanism , Q E3, Q E3 The 

see-saw mechanism 

jisl ] links the low-energy parameters (experimental observables) to high-energy parameters 
(masses of SM singlet neutrinos) through 15 independent parameters of Yukawa matrix. 
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Whichever the direction of prediction by see-saw mechanism goes, top-down or bottom- 
up, the 15 parameters make full contribution in general to the prediction of any physical 
parameters, angles, CP phases, or mass eigenvalues at low energy or at high energy. So far, on 
the other hands, the lack in understanding of Yukawa parameters has restricted consideration 
of the low-high link so as to rely on choice of a model. For instance, there are models which 
predict low-energy observables in terms of parameters of leptogenesis 
and/or minimize the degree of freedom in parameter space of leptogenesis [23 1 
present. 

The purpose of this paper is to find the structure of matrix of neutrino Yukawa couplings 
(Yukawa matrix) which satisfy three constraints simultaneously: (i) the out-of-equilibrium 
condition of leptogenesis, (ii) the upper bound of lepton flavor violating decay in supersym- 
metric Standard Model (MSSM), (iii) the prediction of large mixing angle in atmospheric 
neutrino oscillation j3] and CHOOZ bound j3| of reactor neutrinos through see-saw mech- 
anism. The structure of Yukawa matrix is understood as assembly of a matrix of eigenvalues, 
left transformation and right transformation. An advantage of the above point of view over 
mass structure, to say, in bi-unitary parametrization, compared with taking a whole mass 
structure in weak basis, is that the coexistence of any large mixing angle and hierarchical 
mass spectrum can be clearly defined. Another advantage is that the application of con- 
straints can be systematic because the constraint (i) is irrelevant to the left transformation 
while the constraint (ii) is irrelevant to the right transformation as we will see. While the 
constraints (i) and (ii) have separate dependence on the right mixing and the left mixing, 
respectively, the constraint (iii) is connected to both left mixing and right mixing matrices. 
A series of phenomenological constraints were applied to the three parts of Yukawa structure 
carefully so as to partition parameter space into an eligible part and a ruled-out part. 

First, the out-of-equilibrium condition of lepton number violating process constrains the 
right-handed transformation. In Section II, an appropriate parametrization for model com- 
parison is introduced so that the size of CP asymmetry for leptogenesis can be written in 
a minimal way. In Section III, the out-of-equilibrium condition is rephrased in an effective 
form to fit out purpose and examined model by model depending on leading contribution 
which gives rise to comparison with the bound of the condition. Each model has its own 
characteristic range of elements of right transformation and the expression for the CP asym- 
metry from the leading contribution determined in each model. 
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In Section IV, the strong bound of the branching ratio of fi — > cy in MSSM is briefly 
reviewed. One can derive the constraint (ii) on Yukawa couplings focused on left mixing an- 
gles from renormalization group equation(RGE) of scalar mass terms. As for mixing angles, 
any pattern of left transformation is not ruled out by the constraint (ii). When the consid- 
eration is accompanied with the comparison between eigenvalues, a particular pattern can 
be ruled out by the constraint. The combination of two transformation matrices satisfied 
by (i) and (ii), with the type of eigenvalues, i.e., nearly degenerate or hierarchical, will be 
tested whether they can predict large mixing angle solution in atmospheric neutrinos and 
CHOOZ bound in reactor neutrinos through see-saw mechanism. We have only seven char- 
acteristic combinations which satisfy the three constraints in Tables HT1 and lllll The number 
of combinations is a consequence after excluding the cases where a suppressed eigenvalues 
can be derived only by fine tuning. There are four such cases shown in Appendix A. 2. From 
Table IVl in Appendix A. 2, we will see that the constraint (ii) is redundant if we first apply 
other two constraints (i) and (iii). 

In Section V, we discuss the direction for the improvement of models and the possible 
aspect of application of further constraints in bi-unitary parametrization. A model will be 
examined as an appealing one for leptogenesis. A summary follows. 

II. BI-UNITARY TRANSFORMATION OF NEUTRINO YUKAWA MATRIX 

The leptogenesis has attracted attention recently because the see-saw mechanism is a plau- 
sible way for providing the scale for a lepton number violating process in the early universe. 
The observed neutrino mass spectrum in experiment is consistent with a see-saw model with 
two scales, A EW electroweak(EW) symmetry breaking scale and Aqut the scale of grand uni- 
fied theory(GUT). In the low-energy effective theory, a lepton number- violating interaction 
results in light neutrino masses by non-zero vacuum expectation value of a Higgs scalar, and 
the mass matrix is diagonalized by 3 mixing angles, 1 Dirac phase and 2 Majorana phases. 
As far as the renormalizable gauge invariant couplings are concerned, the lepton sector will 
be 

C = -y E Hle R - y N Hlv R - l -M R W R v R + h.c, (10) 
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in the weak basis. The physical basis of light neutrinos and heavy neutrinos where the 6x6 
mass matrix is diagonal can be obtained by a 6 x 6 unitary transformation from the weak 
basis of left-handed neutrinos in the I and right neutrinos Vr. This full transformation is 
reduced to the MNS j^J in the low-energy limit, in the basis where charged lepton mass 
matrix and right neutrino mass matrix Mr are diagonal. 

Considering that the Yukawa matrix is not necessarily diagonal in the physical basis where 
masses of Majorana neutrinos are diagonal, the transformation of Yukawa matrix (neutrino 
Dirac mass matrix) needs to be defined in addition to the transformation of Majorana mass. 
However, there might exist constraints that can connect the MNS to left-handed mixing 
matrix of Yukawa matrix, for example, arising from a flavor symmetry. Choose a basis 
where M R = Diag (Mi, M 2 , M 3 ). Let 3^d = Diag (yi,y 2 , 2/3), so that a 3 x 3 matrix y N for 
three-generation neutrinos can be expressed in use of bi-unitary transformation as 

y N = hy D Rl (11) 



Pascoli, Petcov, and Rodejohann also discussed this parametrization in Ref.[28j in compar- 
ison with other parameterizations to describe leptogenesis. The unitarity of L implies that 
any left-handed mixing angle in Yukawa matrix does not affect any parameters in leptoge- 
nesis, because the loop contribution appears in terms of y j N y N = Ry 2 D Rl 



ii = y 2 1 \R 11 \ 2 + y 2 2 \Ru\ 2 + yi\Rn\ 2 , (12) 
yfRuRt! + y 2 2 Ri2R* k2 + vlRisKs, (13) 



Ik 

where Rij is an element of R and k = 2 to 3. Expressing the CP asymmetry in Eq.fjHJ) in 
terms of the Yukawa matrices yields, 

2 / , \ 2 



ei « KTTm 



y * yN ) 12 Mi (y^y)^ Mi 
S>t v\ Mo + /S>t^ \ M 



y 1 N y N ) n M2 (ator) 



(14) 



and all the possible models can then be classified into three cases depending on the 
dominant term in Eq. (jl2j) . Factors will be indicated in bold strokes in the following if they 
are much smaller than order one. These lead to the following cases. However, the terms 
with the bold factors cannot be necessarily neglected in the leptogenesis analysis. 
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Case I, where y1\Ri^\ 2 is dominant, 



J/3 1 #13 1 + 2/3 1 #13 1 + 2/3 1 #13 1 

yi|Rll| yii2j ie *» + ^^\y 2 Rl 2 e^ + ysRl 3 e l A , (15) 



,ya|R-i3| ys|Ri3 

where 71,72, and 73 are arguments of Rn,Ri 2 , and i?i 3 , respectively. 
Case II, where j/f I-R12I 2 is dominant, 



yl,y N ' 



,,. _ f y x R n Rl x t y 2 Ri 2 R* k2 y 3 Ri 3 R* k3 



2/2 #12 J/2I-R12I 2/2 # 
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y,|R,l| yi ^ ie ^ + 2/ 2 #* fe2 e i72 + N^ya^a^ , (16) 



,y2|R-i2| " y2|R-i2 

Case III, where yf\Ru\ 2 is dominant, 



,,. _ 1 y x R lx R* kx t y 2 Ri 2 R* k2 y 3 Ri 3 R* k3 



yi\R n \ yi\Ru\ yi\R 



yiR * kl ^ + y^\y 2Rl2e ^ + ZilM^^ , ( i 7 ) 

yi l-rf-ii I yil-ttiil / 



and Case IV, where yf\Rn\ 2 , y%\Ri 2 \ 2 , and y 3 \R X3 \ 2 are equally dominant, 

i y » yN )ik „ fvlRuRll , J/2^12^2 , J/3^13^ 



J/1I-R11I J/2I-R12I J/3I-R13I 
(yii? ie *» + j/ 2 ^ 2 e^ + j/3^ 73 ) 2 • (18) 



There are also cases with two equally dominant terms in ^Ir^jvJ m Eq. ffT^j) . 
Case II*, where y 2 \R 13 \ 2 ,y 2 \R 12 \ 2 > y 2 |#u |2 



ill > 



Case III*, where yf|#i 3 | 2 , y\ |#u| 2 > 2/||-Ri 2 | 2 , 

and Case IV*, where yf\R n \ 2 , y 2 \R 12 \ 2 » yf|-Ri 3 | 2 . 
The corresponding expression to (y^y^ / (y^y^ can be found with a single sup- 
pressed factor. Such 3 possible cases have similar aspect to the cases II, III, and IV when 
the out-of-equilibrium condition of lepton number violating process constrains the magni- 
tude of R, to be discussed further in the next section. 



The hermitian operator 3^y3^v can be expressed with three independent CP phases and 
six real parameters such as 



Zll 



zi 2 e 

-12' Z 2 2 -2:',< 

■y„„,D-i<t>3 



i(f>3 



(19) 



/ 



z i3 e z 23 e z 33 

where the Zij S £1X6 real. Two phases, say 0i and 02, can be eliminated by a diagonal phase 
transformation, P = Diag (1, exp(—i<pi), exp(— ifa))- 

( 



Zl3 



\ 



V 



Z23e 

Z33 



(20) 



/ 



Zll Z12 
z \2 Z22 

Z13 z 23 e~ 1 

where 5' = — 0i — 02 + 03- In general, the diagonalization of Fy^ N y N F^ still requires 3 
independent phases as combination of 2 Majorana phases and 1 Dirac phase. One can take 
the absolute values of the elements of R in Eqs. (fT5j) - (fT7j) to find \yN yMik: regardless of the 
complexity with 71,72, and 73 in Eqs. (|15J) - (fT7j) in terms of elements of R. For example, the 
size of CP asymmetry in Eq. (|14jl can be written as, 
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-MyNW Mi 



y Ny N11 



smzcpi 



Mo 



, n -i \y N y N W . , Mi 

+ 10 — — sm20 2 — , 

M3 



y N yN U 



(21) 



using Eq.(Hnj). 



III. OUT-OF-EQUILIBRIUM CONDITION AND ITS CONSTRAINTS ON 
YUKAWA MIXING ANGLES 



The decay width of Vri by the Yukawa interaction at tree level is 

87T V J 11 



(22) 



The out of equilibrium condition results when the Hubble parameter exceeds the decay rate, 
expressing Hubble parameter in terms of temperature T, 



Ti<H= l.mg l J 2 — 



P l 



(23) 



At temperature T = Mi, the condition can be rephrased as, 

< Ci 2 , (24) 
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where 



9 Mi 

c ' 510 ^- (25) 

The upper bound Ci in Eq.(J2U) cannot exceed Mi/Mqut, when Mqut is chosen in a range 
(10~ 4 — 10~ 3 ) Mpi depending on a theoretical framework and the masses of heavy Majorana 
neutrinos cannot be higher than Mqut- As long as Mi is hierarchically smaller than the 
other masses of heavy neutrinos, as assumed here, the (f should be considered safely as 
small as order of < 10 -2 . 

The condition was examined in detail by Buchmuller and Plumacher 
of generated B — L asymmetry to the effective mass defined as 



121 with the relation 



2 



''' ' l^) llM ? (26) 



Besides the upper bound which keeps the asymmetry protected from wash-out effect, there is 
also a lower bound for mj since, at high temperature, weak Yukawa couplings cannot produce 
enough neutrinos. The eligible fh\ capable of generating sufficient B — L (or equivalent ly 
L) asymmetry, lies in the range mi < rhi < due to see-saw mechanism ^iJIsjL where 
mi is the mass of the lightest neutrino, which is presumed to be non-zero but smaller than 
m 2 ~ a/ Am| and m 3 is the mass of the heaviest light neutrino whose magnitude is about 
a/ Am 2 tm , assuming a hierarchical light neutrino spectrum. 

Using Eq.flHJ), the condition in Eq.(j2U) can be rephrased in such a convenient way as 

y?|i2n| a + IfeW + ySl«i3r<Cf- (27) 
The first significant implication of the above condition is 

l*is|<Ci, (28) 

when y 3 is of order one. The hierarchy in masses of right-handed neutrinos is a preliminary 
for the CP asymmetry in Eq. (JJJ) of leptogenesis, causing the smallness of d and so a small- 
angle constraint for |i?i3|. The next implication is the smallness of yi, 

yt < Ci, (29) 
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TABLE I: Model classification: Cases I, II, III, and IV are introduced in Eas. (jl5H18() . Case -a 

is for y 3 > y 2 > yi, Case -b for y 3 ~ y 2 > yi, and Case -c for vy 3 » y 2 ~ j/i. In |K|, Ci is 
denned in Ea. (|24|) . £ s 's in Cases I-b, IV-b and in Case Il-b represent small values of |-Ri 3 | and 
I-R12I, respectively, while ( s in Case III and Case IV-a represent simply a small angle. The S12 
and S23 imply that those angles can be large or small. The Case II and the Case -c cannot be 
compatible with each other. 



Models 



|K|, 



ei in terms of leading contribution 



I-a,-c 



J/3 2 |«13| 2 »% 2 |«12| 2 ,S/fl«ll| 2 
2/3 > J/2 > J/i or y 3 > y 2 ~ J/i 



1 s 12 < Ci 
512*23 $23 1 



10- 1 fe|R23l) 2 i| 



/ S/3 2 |Rl3| 2 »?/2 2 |«12| 2 -3/ 2 |«ll| 2 
I-b { 

J/3 ~ J/2 > J/1 



1 < Cs < Cl 
< Cs S 23 1 



10- 1 (j/3|i?23|) 2 ^ 



II-a 



j/l|i?i2| 2 »jyi|Ri 3 | 2 ,j/?|flii| 5 

J/3 > J/2 > J/1 



1 S12 < Cl 
S12 1 S23 
5l2«23 «23 1 



10-Hj/2|R22| + *|j/3|R23|} 2 ^ 



II-b 



J/i|Rl2| 2 »J/l|Rl3| 2 ,J/?|Rll| 5 
J/3 ~ J/2 > J/1 



1 <Cl <Cs 

< Cl 1 s 23 

< Cl »23 1 



io- 1 ^!^!) 2 ^ 



J/ 2 |«ll| 2 »J/ll«12| 2 ,J/i|fll3| 2 
III-a,-b,-C { J/3 > J/2 > J/1 , J/3 ~ J/2 > J/1 , 
or J/3 > J/2 > J/i 



1 G <Ci 

Cs 1 s 23 
Cs S 2 3 1 



+ ^& 2 |R22| + ^>|i?23|} 2 ^ 



IV-a,-b 



y 2 3 \Ri 3 \ 2 ~ y 2 2 \Ri2\ 2 ~ yf\Rn\ 2 

J/3 > J/2 > J/1 or J/3 ~ J/2 > J/1 



1 C s <Ci 

Cs 1 S23 
Cs S23 1 



lO-itel^l+j/sl^sl} 2 ^ 



IV-c 



J/l|Ri3| 2 ~jy 2 2 |Hi 2 | 2 ~j/ 2 |/?ii| 2 

J/3 > J/2 ~ J/1 



1 1 <Cl 
1 1 S 2 3 
S23 «23 1 



10- 1 {j/2|R22|+J/3|iJ23|} 2 ^ 



since \Rn \ is always of order one whether mixing angles are large or small. So y\ cannot be 
of the same approximate size as 2/3 if 2/3 ~ 1. The models that satisfy the two conditions 
in Eqs. fl2%j) - (J2§j) can be according to which term dominates in the real number (34r3^v)n 
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and how different the eigenvalues yi,y2, and 2/3 of the Yukawa matrix are. A model with 
2/3 y2 ^> 2/i is denoted as Case -a, the one with 2/3 ~ 2/2 ^> y\ as Case -b, and the one 
with 2/3 ^> 2/2 ~ 2/i as Case -c. The leading contribution to e\ in Table is presented by the 
first term in Eq. (jl4j) . 

In Case I in Eq.(|15|). if the eigenvalues are of normal hierarchy, y 3 ^> y 2 , the out-of- 
equilibrium condition in Eq. ()27|) does not place any constraints on elements in M other than 
I-R13 1 < Cij while nearly degenerate eigenvalues, 2/3 ~ 2/2, restrict | -R12 1 to be small. If the 
transformation M./JJ = MJ _1 is parameterized as 



1 C13C12 S12C13 S13 

~ s 12 c 23 — S23 s 13 c 12e lS C23C12 — S23 s 13 s 12 e lS S23 c 13 e tS 
y S23Si2e l5 — S13C23C12 — S23Ci2e lS — S13S12C23 C23C13 J 



(30) 



where J is a diagonal phase transformation, one can deduce the smallness of I-R21I and I-R31I 
from the smallness of Ru due to small S\2- The possible structures of K for the Case I are 
listed in Table HI where S12 and S23 specified in the table are allowed to be either small angle 
or large one, while X s indicates that only a small angle is allowed. If S23 is not so small, the 
leading contribution of the imaginary part which leads to e\ consists of mainly (2/3I-R23I) 2 
and M1/M2. The model examined in Ref . \v\ is an example of Case I-a. 

In Case II, in Eq. (jl6|) . if the eigenvalues are of normal hierarchy with y 3 ^> y 2 , the 
out-of-equilibrium condition allows I-R12I to be much larger than \Ris\. If the assumption 
is made that 2/2 is not larger than (1, then I-R12I can still be a large angle. In other words, 
I-R12I can be either small or large, though, it is larger than (j in this model. The leading 
contribution to e x depends on mainly from 

y*3 1 j^,^ 2 1 

2/2I-R22I H P5 — r 2/3 1 J?23 1 - (31) 

y2|-tVL2| 

Otherwise, the two eigenvalues, y 2 and 2/3 in Case II, are close to each other. The element 
I-R12I is smaller than £1 and |i?i3| is far smaller than \Ru\- The term |3^jv3^Af|i2 / (y^y^ 
reduces to 2/2I-R22I which is now the leading contribution to e\ of order one. The assumption 
implies I-R12I ^> 2/i/2/2> so that the Case II-c with 2/1 ~ 2/2 is ruled out. For the same 
reason, there can not be Case II*-c, while Case II*-a with y 3 S> 2/2 and Case II*-b with 
2/3 ~ 2/2 and small I-R12I are allowed. The allowed ranges in the elements of |R| for the Case 
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II*-a and Case II*-b can be obtained as the same as in the Case Il-a and Case Il-b, 

respectively. 

Case III in Eq.(JT7J) corresponds to |i?i 2 | <C 2/1/2/2, i-e., I-R12I is always a small angle. Any 
part of the contribution to |3^jviVjv|i2 / ^3^Jv3^j\rJ is not yet ruled out by the assumptions 
given here. The Case III* also gives rise to similar constraint on I-R12I as Case III does. 

In Case IV in Eq.()18|). a subcase with j/3 ~ y 2 and I-R12I is not allowed by out-of- 
equilibrium condition because 2/2I-R12I should not be larger than £i- This is the case with 
~ 2/1/2/2- Thus, |i?i 2 | is small for 2/2 ^ 2/1, while |i?i 2 | is large for y 2 ~ 2/1 • The possible 
structures of |M| derived in Case IV* also is consistent withe those in Case IV. 

The structures of |R| in Cases II*, III*, IV* attribute their correspondences to those 
in Cases II, III, IV, respectively, to relative significance in comparison of 2/1I-R11I and 
2/2I-R12I when y 3 \Ri 3 \ is fixed. In other word, once I-R13I is bounded, cases can be effectively 
classified into a case with 2/1I-R11I 2/2I-R12I, a case with y 2 \Ri 2 \ "C yi\R\±\, and a case 
with ~ 2/2 1 -R12 1 • Thus, the cases with two equal dominant terms can remain without 

further detailed discussion. 

In Table HJ where the various cases are considered, the allowed structures of |M| under 
those cases and the CP asymmetry ei in terms of leading contribution after ruling out the 
suppressed part under the assumptions are summarized. In Case I and Case IV with I-R23I 
of order one , or in Case Il-b, e x can reach its model independent maximum value, 

-i M i 

Clmaz ~ 10 — , (32) 
1V12 

where M 1 < M 2 . 



IV. LOW-ENERGY OBSERVABLES AND THEIR CONSTRAINTS ON 
YUKAWA MIXING TRANSFORMATIONS 

The current experimental limit on the branching ratio {BR) of the lepton flavor violat- 
ing(LFV) decay mode, \i — > is 

BR -> erf) < 1.2 x 10 -11 . (33) 

In the SM, the branching ratio is < 10~ 50 suppressed far below the observable bound |29|. 
Another theoretical framework where one can estimate the size of the LFV decay rates 
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is the minimal Supersymmetry(SUSY) Standard Model(MSSM), where the rates for LFV 
processes can be enhanced due to large tan/? [s(3] . The neutrino Yukawa couplings cause 
the renormalization group equations (RGE) to develop flavor changing contribution to soft 
masses of scalar leptons, which may originate at a universal value for all kinds scalar leptons 
at GUT scale Mx- The generated LFV mass terms for scalar leptons after integrating RGE 
is 

(A-!).- - ^(3^ + ^)£3Wog§^ 

k k 

-1 o .ox. My 



(Zml + Al) log ^(y^)^ (34) 



8ir 2 * " ivi 3 \ " / ij 

where log j^- for k = 1 — 3 can be considered as of the same order even when M^'s are 
of hierarchy. The tjIq and A$ are the universal masses and universal trilinear couplings, 
respectively, in soft SUSY-breaking lagrangian 2l3|. The flavor changing decay in Eq. (j33j) 
involves neutralino exchange and chargino exchange in loop diagrams, and the resulting 
amplitude is proportional to tan/3. Therefore the branching ratios can be expressed in 
terms of 3^v3^iv as follows; 



BR^^l^)- (35) 
" 3 ' _1 (3m^ + ^)log^(^)j 2 tan 2 /3, 



Glm* 1 8tt 2 



where m s is the typica 



mass of a superparticle, a is the fine structure constant, and Gp 
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3ll. Recently, LFV in low-energy has been considered for its 



is the Fermi constant 
relevance with leptogenesis 

Even if a model can predict an enhanced branching ratio near the range accessible in a 
near future experiment, the eligible structure of Yukawa matrix implies there might be a 
constraint for severely suppressed mixing angles in left-handed transformation of Yukawa 
matrix at low energy. The enhancement mechanism of branching ratios is possible not only 
by raising tan/3 but also by increasing mixing angles. Hisano et al H- analyzed the CKM 
matrix of quark mixing. With small mixing angles as in CKM, the branching ratios of 
— * £jl processes for large tan/3 can reach close to the current experimental bounds in 
Eq. (}3*3*j) . Considering the strong upper bound in current experiment limit for the decay 
mode fi — > erf, the relevant constraint on the Yukawa couplings is, in terms of the elements 
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of L defined in En.(fTT|). 



y N y 



N 
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yjL 2 iL* n + y\L 22 L\ 2 + ylL 23 L* 13 < 1, 



(36) 



which can be satisfied only if y 2 or |L 12 | is small and |L 2 3| or |L 13 | is small. 

In the rest of this section, I combine the constraint on R from the out-of-equilibrium 
condition and the constraint on L from the upper bound of the branching ratio of p, — > cy. 
The question is whether the result of those combinations can be accommodated with the 
experimentally observed large mixing angles of light neutrinos when implemented via see- 
saw mechanism. Light neutrino masses can be obtained through the see-saw mechanism 
where Dirac masses are defined in terms of the non-zero vacuum expectation value v 2 of a 
light Higgs in MSSM with two Higgs doublets, 



m v = -v 2 2 y N M^yl } 



lhy D R^M^R*y^h T . 



(37) 



In this basis the mass matrix of charged leptons is diagonal. When the transformations L 
and M are trivial, i.e., equal to the identity, the see-saw mechanism yields the light neutrino 
mass matrix, 



vl vl vl 



-v 2 Diag — , — , — 



(38) 



Mi M 2 M 3/ 

from Eq. (j3*Tj) . It will be interesting to watch how neutrino mass matrix in Eq. (j37j) can change 
as the type of transformation matrix M, and subsequently L, switches from one to another 
general form. There are 4 types of K's, which have been given in Table HI and 4 types of L's 
are as follows: 



|Li| 



/ 



1 A A 
A 1 1 
A 1 1 



\ 



1 1 p 



|L 2 



1 1 A 
1 1 A 
v AAly 



(39) 
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which satisfy the out-of-equilibrium condition in Eq. (J24|) . The most simplified notation is 
used, which is that diagonal elements and large mixing elements are denoted by 1, whereas 
small mixing elements are denoted by arbitrary small values by p <C 1 and A <C 1 in R. and 
L, respectively. That is, there are only two kinds of elements in a transformation matrix: 
order one or far less than order one. The more detailed Ws have been sorted in Table U 
by considering out-of-equilibrium condition. Later, in case the leading order in an entry 
of light neutrino mass matrix originates from small angles in transformation matrices, the 
light neutrino mass matrices in terms of the small angles, considering the actual entries as 
expansion of small angles, are shown for models in Appendix A.l. 

Each structure of |Lj| — |L 4 | is similar to each of |Ri| — as shown in Eq. (JH§|) . which 
satisfy the condition in Eq. (j36|) which is the same form as Eq. (j27|) for the M. matrix. For the 
suppressed (3^3^)21, it is necessary that L 12 and L 13 L 2 3 be small for y 2 ~ 1/3, while the 
smallness of L 12 or L 23 is optional for y 2 ^ and L 13 small. The possibility of large |L 13 | 
is excluded by the CHOOZ bound [26| on the MNS neutrino mixing matrix. 

In Tables |H] and IIII[ all the possible structures of light neutrino mass matrix, which can 
be produced by Eq.^ using Ws and L's in Eq.flHHJ), are listed. Due to M x < M 2 ,M 3 
assumed for Eq.((7j), the following approximation can be used in the tables. 

Mf 23 ~ Mf/ ~ Mf 1 , 

Mr 2 1 3 »M 23 1 , M 1 - 1 3»M 3 - 1 , (40) 

where 

M^ 2 \ = M^ l + M 2 l + M^\ 

M£ = M? + M;\ (41) 
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TABLE II: List of possible structures of matrix m v jv\ which can be produced by see-saw mechanism 
in Ea. ()37|) using R's and L's in Ea. (|39|) . where m v is a matrix mass of light neutrinos and V2 is a 
vacuum expectation value in two-Higgs-doublet models. M123 and M12 are defined in Ea.(|41|l and 
yi23 and yx2 are defined in Eq.([42|). Please see Appendix A.l for the entries marked by *. 



Li 



La 



yjM~ 



y 2 23 M 23 vh M 2s 



y 2 l M 1 1 +y 2 2 M 23 1 y 2 M 1 1 +y 2 M 23 1 y 2 y 3 M 23 x 
2/22/3M" 1 Mr 1 M" 1 



•|M. 
y 2 y 3 M 23 



y\M 23 



vfM- 1 yiy 2 M- 2 1 yiy 2 M^ 
V1V2M- 1 ylM-^+ylM- 1 ylM-^+ylM- 1 
yx V 2M- 2 x 2/f M -i + y 2 M -i y 2 M -i +3/ 2 M -i 



y\2^i2 

y\2 M 12 y\2 M 12 



Hi 



yl^ 2 \ 



yiV2M 12 \ + yiy 3 M 23 yiy 2 M 12 \ + yiy 3 M 23 
yiV2M- 2 \ + yryaM^ 1 y\M^ 2 \ + y 23 y 3 M 23 x y\M~ 2 \ + yizyzM^ 

r—l 1 tf-1 ..2u-l 1 tf-1 ..»»<■- 1 1 .. .. 



3/i2/2M 12 3 + yrysM^ 1 y\M X2 \ + y 2 3y 3 M 23 1 



vt^+y^M^ 



vi 2 M i 2 \ y\ 2 M x 2 \ yi2y 3 M 23 

y 2 i 2 M i 2 \ y 2 u M i2\ 
y\2yzM 23 y 12 y 3 M 23 y 3 M 23 



• 2 ,M-,„\ y\ 2 yzM 23 

1 rlM" 1 



y\M^ 



y\M- x + ylM- 1 ylM- 1 + y^M^ 1 
Vl M 2 1 + y\ M 3 ~ 1 y'i M 2 1 + y'l M 3 1 



1 +y 2 2 M 2 1 y'lM^+ylM- 1 
y\M^ x + y\M 2 x y\M~ x + y%M 2 x 



y\M 3 



"Model A 
ft Model B 
c Model C 



(Ri,Li), y 2 -V3 
(M 2 ,Li), y 2 « y 3 

(K4,Li), y 2 - Vs Model D: (R^Li), y 2 < Vs 



etc., and 



2/123 = Vi + 2/2 + 2/3, 

2/12 = 2/1 + 2/2, 2/23 =2/2+2/3- 



(42) 



The comparison between M 2 and M 3 is not yet fixed. Either possibility, M 2 <C M 3 or 
M 2 ~ M 3 , is still considered. 

The structures listed in Tables |H] and IIHI are investigated whether they can predict two 
large mixing angles and one small mixing angle in light neutrinos. As for the large mixing 
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TABLE III: The matrix m v jv\ which can be produced by see-saw mechanism in Ea.(|37|) using R's 
and L's in Ea. (|39|) . The % in the matrix from R4 and L3 runs 1 to 3. Please see Appendix A.l for 
the entries marked by *. 



L 4 




yfM- 1 



V2V3M- 1 yjM- 1 



y\i M vi 



i-llz 



y\i M xi y\i M ™ +al M3- 1 y\ 2 M 12 x 



■vi M 3 



ylM- 1 



y\2 M \?i y\2 M \?A + J/wifeM^ 1 y\ 2 M^ 2 \ + yi^M^ 1 

y 2 i2 M u3 + 2/122/3M- 1 y\ 2 M- x + yj^M- 1 yl 2 M~ x + jft m M£ 
yi2 M i23 + yizyaM^ 1 y\ 2 M- x + y 2 12S M 2 ^ y\ 2 M~ x + y\ 2i M 2Z x 



2/? M l~23 VlV2 M vk 2/1 2/3 M" 1 

2/12/2 iWi" 2 3 y\M~ 2 \ 2/22/3 M^ 1 
2/1 J/3 M 2 3 WzV3M^ y\M 2Z x 



y\M- x + y\M 2 x y\M~ l + y\M 2 x yjAff 1 
2/ fM 1 - 1 +2/lAf 2 - 1 j/fMr 1 



y\ Mf 



2/1 M^ 1 



2/1 



"Model E: (Ki,m),J«!~I/3 
^Model F: (R 2 ,L 3 ), 2/2 < y 3 
c Model G: QR 4 ,L 3 ), y 2 « y 3 



angle in atmospheric neutrino oscillation, the ratio of two elements m^/m^ must not be 
less than order one. On the other hand, for the ratio m v \ijm v ^i to be not less than order one 
is a sufficient condition for the large mixing angle in solar neutrino oscillation, because there 
is still a room for a kind of manipulation including fine tuning so as to make the mixing angle 
large even without the naive ratio m u i2/mv22- The application of the constraints from light 
neutrino mass spectrum is described more in Appendix A. 2. When y 2 ~ 2/3, the combination 
of R and L which help us predict large mixing angle for atmospheric neutrinos is (Mi, Li) for 
Model A, (R 4 , Li) for Model C, or (R h L 4 ) for Model E where R { and Lj with i = 1 - 4 
are denned in Eqs. (j39|) . The solutions R2 or R 3 are not included since the out-of-equilibrium 
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TABLE IV: List of phenomenologically viable models of Yukawa matrices and the light neutrino 
mass matrix divided by v 2 derived through see-saw mechanism. The eigenvalues in Models A, 
C, and E are of y2 <C 2/3 ~ 1, while those in Models B, D, F, and G are of y\ <C j/2 ~ 273 ~ 1- 
The 2/2 and yz specified in m u /v 2 remain for a later use. 

Model |L| y D |Rt| m„/v 2 ~hy D R J <M R 1 M.*y'[,l, T 
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A 




yi 
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p 
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1 




1 
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1 
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1 




1 




_ p 
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1 





yfM- 1 v v 

(ylM- 1 +yl. i M^)\ + y 1 y 2 ?,M w \p y^M^ 1 yj 
(yjM- 1 +y% 3 M- 1 )\ + y 1 y 23 M- 1 3 p y^M' 1 y^M" 1 
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p 
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1 
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ylM- 1 ^ V 

yxyiM^ y^M^ 1 +y? i M.- 1 y/ 
yxyiM^ y\M-^ +y? i M.- 1 y^M^ 1 + y^M^ 1 



1 A A 
A 1 1 
A 1 1 



yi 



1 p p 
pip 
p p 1 



ylM- 1 v V 

y 2 l M7 1 \ + yi ( y2+yi )M- 3 1 p ylM^ 1 + ylM^ 1 y/ 
y\M7 x \ + yi (y 2 + y 3 ) M^p ylM' 1 + y'iM- 1 ylM' 1 + y^M' 1 



D 



1 A A 
A 1 1 
A 1 1 



3/1 



3/2 



1 p p 
p 1 1 
p 1 1 



y\ Mr 1 yj V 

y 2 M7 1 \ + yi (y 2+ y 3 )M- 3 1 p ylM^ 1 + ylM^ 1 ^ 
y\M- x \ + j/i (y 2 + y 3 ) M'^p y^M' 1 + y\ M' 1 j^M" 1 + yjM' 1 
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p 
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yfM- 1 

(ylM- 1 + y 2 y 23 M- 3 l )\ + y^M^p y\M^ 
(y\ Mf 1 + j/3 3/23 JWm 1 )A + yiya Mf 23 P 3/2 3/3 M^ 1 



3/2 3/3 M^ 1 3/fM^ 1 
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1 


A 




3/1 
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P 
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1 




3/2 
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P 
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1 




1 
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1 





y\ 2 M^ y^M^+ylM.- 1 V 

y\i M x 2 y 2 i2 M i 2 +yi M s 1 yh M i 2 + y 2 ^ 1 



G 



1 1 A 
111 
111 



3/1 



3/2 



1 p p 
pip 
p p 1 



ylM- 1 + y^M- 1 J J 
ylM- 1 + ylM- 1 yjM' 1 J 
ylM- 1 + ylM- 1 yfMr 1 y^Mr 1 
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condition with y 2 ~ 2/3 ruled out a large angle at R12, as shown in Table |TJ When 2/2^2/3, 
(R 2 ,Li) for Model B, (R 4 M) for Model D, (M 2 ,L 3 ) for Model F, (R 4 ,L 3 ) for Model G 
can predict the large mixing angle for atmospheric neutrinos and CHOOZ bound of reactor 
neutrinos. See Appendix. 

Thus, the structures of Yukawa matrices that satisfy three constraints: 

(i) the out-of-equilibrium condition of lepton number violating processes, 

(ii) the upper bound of the branching ratio of the rare decay /i — > ej, 

(hi) the large mixing angle in atmospheric neutrino oscillation, 

must be one of the 7 combinations, Models A to G in Tables ITT1 and lllll which establish the 
relative dominance of y% and 1/3, and the approximate magnitudes of transformation matrices, 
|R| and |L|. In Table ITVl is the list of models of Yukawa matrices and the light neutrino 
mass matrix divided by v 2 . Models are specified in terms of two unitary transformations 
and the relative dominance in eigenvalues as summary of Tables |H] and IIHI The notations 
are introduced in Eqs.(|39~ |) .pT ]) .pS ]) . and (jAlj) - (jA3j) . The entry marked with y/ represents 
the corresponding symmetric elements. There are two types of spectrum of eigenvalues, 
Z/i <C 2/2, yi ~ Jj2 *C JJ3 ~ 1 and y\ <C y% ~ 2/3 ~ 1 , which correspond to Case -a,-c and 
Case -b introduced for Table [IJ respectively. 

In the Appendix, Table IS describes in detail which constraints keep certain entries in 
Tables ITD and Hill from being eligible models. The Table IVl in Appendix A. 2 shows obviously 
that the cases prohibited by the constraint (ii) are already prohibited by the constraint (iii). 
The constraint (ii) turns out as a redundant one in checking the eligibility of Yukawa matrix 
to known leptonic phenomenology. In other words, the structure of Yukawa matrix which 
can derive the light neutrino mass matrix through see-saw mechanism always appear to give 
rise to Br(fi — ► ej) below the present experimental limit in SUSY theories with large tan/3 
and universal slepton mass at GUT scale. It is possible to pull out a constraint to narrow 
down the estimation of heavy neutrino masses from LFV through the see-saw mechanism. 
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V. MODELS AND A LINK BETWEEN LEPTOGENESIS AND LOW-ENERGY 
NEUTRINOS 



A. Application of bi-unitary parametrization 

The light neutrino mass matrix defined through the seesaw mechanism in Eq.(|37j) is diago- 
nalized as 

Vm L V T = -v^Ly D R^ M^-Wy^L T , (43) 

where V is MNS light neutrino transformation matrix and is a diagonal matrix with 
mass eigenvalues m^. The equation of 3 x 3 symmetric complex matrices in Eq. (j43|) consists 
of 12 equations of parameters, while its right hand side consists of 18 unknowns: 3 angles 
and 3 phases in L, 3 angles and 3 phases in R, 3 eigenvalues in y^, and 3 heavy neutrino 
masses in Mr. Once one makes a choice out of the 7 models model A to model G specified 
in Tables |H] and IIII1 one can draw the boundary of eligible parameter space starting the 6 
angles, large or small. Even though it is supposed that the number of the free parameters 
reduces to 12 for 12 equations, considering 6 angles in |R| and |L| fixed, solving Eq. ()43)1 
yet imbeds a great deal of ambiguity, since CP phases in MNS matrix and individual mass 
eigenvalues in left-hand side can barely be said as fixed. 

There are a few approaches to make improvement in solving the see-saw mechanism 
depending on models. First, one can rely on theories with additional symmetries to reduce 
the number of parameters. For example, in case of Model G, a flavor symmetry like 
U(l) can be utilized to simplify the mechanism in Eq. ()43j) . because the charges of the 
flavor symmetry assigned to vm (£i) can result in the common constraint on the side of 
right (left) mixing in Yukawa matrix and the structure of right (left) neutrino mass matrix. 
The structure of Yukawa matrix in Model E can fit in a theory of grand-unified-type 
symmetry between leptons and quarks, if Li is close to CKM matrix. Among eigenvalues 
of Yukawa matrix and 3 heavy masses, there are a few parameters which can be fixed by 
hands motivated from a higher-rank symmetry. They can be of order one and/or M3 of 
GUT scale. The 2/3 ~ 1 has been used in a number of times since Eq. (j27j) . as well as 1/2 ~ 1/3 
in particular models. 

Another possible approach can be collecting more phenomenological constraints involved 
with neutrino Yukawa couplings to diminish the number of the eligible structure of Yukawa 
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matrix, aiming the survival of only one natural structure. The classification derived based 
on bi-unitary transformations is useful in both approaches. That is, the models parame- 
terized using bi-unitary transformation of Yukawa matrix can be easily tested for the com- 
patibility to additional symmetries as proposed in the previous paragraph. The bi-unitary 
parametrization is useful also to examine the eligibility to additional phenomenology. In 
describing a process, the loop contributions of Yukawa couplings consist of a hermitian op- 
erator y^y or yy*. It is always possible to rephrase operator in terms of the minimal 
number of parameters, the pair of (^d>K) or the pair of (3^3, L) without any mixing of 
the two transformations. Thus, application of an additional constraint will be simply the 
re-examination of the 7 models characterized by eigenvalues of Yukawa matrix, left and right 
transformations. 

B. A model of neutrino masses 

The description of leptogenesis requires fine understanding of CP phases in Yukawa matrix 
as well as CP phases in light neutrino mass matrix in order to secure a sufficient amount 
of CP asymmetry from the decays of heavy Majorana neutrinos. In lepton sector, however, 
it is not practical to regard the CP phases at low-energy as known parameters. Rather, 
in analogy of quark sector, the size of a real mixing angle can be realized as the source of 
certain amount of CP violation before the derivation of detailed imaginary parts. Many 
models proposed so far took simply order one contribution from CP phase. 

As far as the detailed discussion on CP phases is set aside, Models A, D, and E with a 
large angle at R23 can provide optimal scenario for the sufficient amount of CP asymmetry, 
accordingly lepton asymmetry, Eq.fjSJ). There are a number of cases listed in Table E which 
can enhance the CP asymmetry e\, if I-R23I becomes close to one, so as to reach its maximum 
value. 

« ~ ur'iSHWf ~ 10-^. (44, 

It is worthy of attention that large mixing angles in low-energy originate entirely from the 
large mixings in heavy neutrinos in Model E. 
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The Yukawa matrix in Model E is 

/im\/ 



|3^| = 



1 A A 
A 1 A 
A A 1 



\ 



J 



p i i 



(45) 



yi o o 

1 
1 

where y 2 ~ y$ ~ 1, 2/1 *C 2/2- It is possible that the Model E reduces to the model 
with two-generation of heavy neutrinos as the large mixing angles become maximal and the 
degeneracy in eigenvalues becomes exact. In such limit, the leading contribution to e\ does 
not any longer come from R23. The matrix in Eq. (j45|) was introduced in the basis the mass 
matrix of righthanded neutrinos is diagonal. Now, a new basis can be chosen in such a way 
that Yukawa matrix is symmetric. Let the Yukawa matrix in Model E in the new basis be 



N\ 



1 A A 
A 1 A 
A A 1 



Vi 
1 
1 



(46) 



Then, the righthanded neutrino mass matrix in Eq. 
mixing angle in the new basis, 



p' 1 ft 

V p' V 

will have a structure with a large 



D* TO>T 



(47) 



where IR ^ is the right transformation to diagonalize the Yukawa matrix from the old basis 
in Eq. (j45|) . while M. new is the one for the diagonalization from the new basis in Eq.([46p. The 
Model E can represent a model in left-right symmetry, in other words, the model shows a 
well-balanced correspondence between low energy and high energy over see-saw mechanism, 
since the Yukawa matrix is symmetric and two Majorana neutrino mass matrices are very 
much of analogy. 

The difference of the lepton sector from the quark sector was recognized from the large 
mixing angles. In this scenario, the dissimilarity in the two sectors can be understood as 
originated from the existence of massive singlet neutrinos, whereas nothing like that is in 
the quark sector. There might have been dynamics of heavy neutrinos in the early universe 
to generate the large mixing angles of which become later the source of the large mixing 
angles of light neutrinos. It is shown in Eq . (J44j) that the leptogenesis is optimized. The 
Eqs. ()46)) -([47 )) help the see-saw mechanism look more meaningful with this model in a theory 
of left-right symmetry. 
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C. Leptogenesis and neutrino large mixing angles 



Before concluding the section, we give a some thought on the role of the see-saw mechanism, 
expressing CP asymmetry in terms of light neutrino parameters in a model-independent 
way. Using solutions to Eq.(|43|) for a chosen model, one can express the CP asymmetry in a 
proper way to replace the heavy neutrino masses in terms of low energy neutrino observables 
and parameters in Yukawa matrix. The diagonal neutrino mass matrix in Eq. (j43|) can be 
rephrased as 



(48) 



where 3^v = 3^vIP* is defined in such a way that phases in P do not appear explicitly in y^ 
as hidden in Eq. (j20|) . Alternatively, 



Mk ¥ k 2 = -viy^ — (\ T y N 



the imaginary parts of whose inverse are obtained as 



k'k 



(49) 



sin 20! 

Mo 



sin . 



M, 




^ T y N 



J k2 



-1 



-1 



^ T y N 



k?> 



(50) 
(51) 



while 



-T 1 



Y T y 



N 



kl 



(52) 



When the value of the rectangular brackets of the above equation is of order one, the term 
with m^ 1 is guaranteed to be a significantly leading contribution compared to the other 
terms with m^ 1 or mj 1 inside the summation, resulting in the enhancement in the scale 
of Mi. In Ref . \v\ . it was examined that large mixing angles in MNS matrix can predict 
the small dilution mass preferred by the out-of-equilibrium condition in a model. In the 
right hand side of Eq. (j21|) . sin20i (M1/M2) and sin 202 (M1/M3) can be re-parameterized 
only in terms of Yukawa matrix and light neutrino masses and mixings using Eqs. ()51|) and 
(I52|) . The enhancement of Mi contributed from large mixing angles in Eq. ()52j) results in the 
enhancement in the amount of the CP violation in Eq. (j21j) . 
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VI. SUMMARY 



The CP asymmetry and the out-of-equilibrium condition of lepton flavor violating process are 
analyzed in terms of eigenvalues and right mixing angles. In the out-of-equilibrium condition, 
there are significant implications. The mixing element |M. 13 1 should be small with an upper 
bound, Eq.pSjl. and two eigenvalues y\ and 2/3 cannot be of similar size, Eq.(J2Hj). Depending 
on which term in expression of the condition actually leads the condition or depending on 
whether eigenvalues are of hierarchy or of similar size, the magnitude of right transformation 
is characterized differently, as shown in Table |U Examining all the possible models, we 
found the model independent maximal CP asymmetry to be 10 _1 M 1 /M 2 , Eq . (|32j) . Besides 
thermal out-of-equilibrium condition, we examined also the upper bound of a LFV decay and 
mixing angles in neutrino oscillations to constrain all kinds of parameters, i.e., left and right 
mixing elements and eigenvalues. We found that only the seven models can characterize 
phenomenologically eligible Yukawa matrices. 

The constraint from the LFV did not exclude any choice of parameter region addition to 
the category excluded by two other constraints. In other words, the survival of only seven 
models is the outcome of the application of two constraints: the out-of-equilibrium condition 
and data of neutrino oscillation. In Table the determination of eligible models from the 
application of certain constraints is specified. The bi-unitary parametrization is useful for 
model test of the compatibility to additional symmetries and the eligibility to additional 
phenomenology. The Model E is appealing in a sense that the size of the CP asymmetry is 
maximized and the number of free parameters can be reduced based on left-right symmetry. 
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APPENDIX A: REMARKS ON TABLES II AND III 

1. Recovery of small angle contribution in transformations 

In Tables |H] and 11111 the elements of matrices determined from small angles in transforma- 
tions are not specified. (They are indicated by a *.) The Models A, C, and E include 
elements, the leading order of which is first order in A or p, while other models have their 
leading order all zeroth in A and p. The small values are not all of the same order of 
magnitude in a transformation matrix. However, different small angles do not distinguish 
because the comparison of their sizes is not necessary to examine the eligibility of a model 
to given constraints in this paper. With p's and A's in M's and L's, respectively, the possible 
symmetric structures of matrix m v /v 2 of Models A, C, and E will be as follow; 



Model A 



Model C 



Model E 



V V 
V 



d 1 — 1 V 

(^Mf 1 + y 2 23 M 23 l ) A + ymsM^lp y 2 23 M 23 x 
(y 2 M^+y 2 23 M 23 1 )\ + y 1 y 23 M u 1 3 p y 2 23 M 23 l y\ 3 M 23 l 



yfMr 1 V v 7 

yfM^X + yi (y 2 + y 3 ) M^p y\M 2 x + y 2 M 3 Y J 
yfM^X + yi (y 2 + y 3 ) M£p y\M 2 l + y 2 M 3 l y 2 M 2 l + y 2 M 3 x 



yfMr 1 V V 

(y 2 M^ + y 2 y 23 M 23 x ) A + ymM^p y\M 23 x yj 
(y 2 M{ 1 + y 3 y 2 3M- 1 ) A + ymM^p y 2 y 3 M 23 l y 2 M 23 



(Al) 



(A2) 



(A3) 



where i = 1—3 and the entry marked by a/ represents the corresponding symmetric elements. 

The 1-2 sector in light neutrino mass needs to be examined to see whether a model is 
consistent with for the large mixing angle of solar neutrino oscillation. In order for a mass 
matrix to have the large mixing angle between the first and the second generations, the 1-2 
element of the matrix m v \ 2 should be not smaller than the 1-1 element m v \\. If one takes a 
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TABLE V: The specification of constraints which in fact rule out certain types of Yukawa matrices 
represented by a 1, a L, and the comparison of y% and 7/3. The CI, C2, C3, and C4 are the indices 
of the constraints explained in Section A. 2. For example, Yukawa matrix characterized by Ri,Li 
and ?/2 *C 2/3 are forbidden due to the constraint C4. 









1L<2 


JL, 3 


If , 


Mi 


2/2 < 2/3 


C4 


C3 


C4 


C3 




2/2 ~ 2/3 


Model A 


C2 C3 


C2 C3 


Model E 


R 2 


2/2 < 2/3 


Model B 


C3 


Model F 


C3 




2/2 ~ 2/3 


CI 


CI C2 C3 


CI C2 


CI C3 


R 3 


2/2 < 2/3 


C4 


C3 


C4 


C3 




2/2 ~ 2/3 


CI 


CI C2 C3 


CI C2 


CI 


R 4 


2/2 < 2/3 


Model C 


C3 


Model G 


C3 




2/2 ~ 2/3 


Model D 


C2 C3 


C2 C3 


C3 



hierarchical neutrino mass spectrum mi ~ a/ Am| and 777.3 ~ V Am\ tm , then the structure 
of neutrino mass matrix should require 



m vV1 / Ami 

777^22 V Am ltm 



"0 



(A4) 



i.e., the ratio between the two elements should be suppressed by one order of magnitude 
|33| . The m v jv\ matrices for the three models in Eqs. (jAl|) - (jA3|) show that there can be 
parameter region which allows large mixing angle solution for solar neutrino oscillation as 
well as the large mixing angles of atmospheric neutrinos and satisfies the CHOOZ bound of 
the reactor neutrinos. 



2. Application of constraints 

The exclusion of certain types of Yukawa structures is explained in Table |3 by specifying 
the constraint violated by each particular structure. The constraint CI is equivalent to the 
constraint (i) in listed in Section IV. The out-of-equilibrium condition of the lepton number 
violating process, Eq.(|24")l. does not allow a structure of Yukawa matrix with both y% ~ y 3 
and R12 ~ 1. 
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The constraint C2 is equivalent to the constraint (ii) in Section IV. The suppression in 
(yNylfj for BR(fi — > ej) < 1.2 x 1CT 11 does not allow the structure of Yukawa matrix 
with both y 2 ~ y 3 and L 12 ~ 1. 

The constraint C3 is equivalent to the constraint (iii) in Section IV. The phenomenological 
constraints for light neutrino experiments require the structure of mass matrix to have 
m vl3 < m u23 ~ m u33 . 

The last constraint C4 is the requirement for Yukawa matrix not to involve fine tuning. 
It is not appropriate to assume any possible fine tuning in Yukawa structure since the 
comparison of parameters was done whether it is hierarchical or nearly degenerate, whereas 
it is unavoidable to assume some possible fine tuning in consideration with light neutrino 
masses due to the coexistence of large mixing angles and hierarchical mass spectrum in its 
symmetric structure. In the Yukawa matrices, if y 2 <C y 3 when L 23 ~ 1 and R 23 ~ 1, such 
a small eigenvalue y 2 can be obtained only by fine tuning. For that reason, the cases with 
L 23 ~ R 23 ~ 1,2/2 ^2/3 are ruled out. In actual model construction with respect to data, 
however, the structures restricted by C4 can be regarded rather positively, if the fine-tuning 
is the only way to accommodate data. If one allows tuning with L 23 ~ R 23 ~ 1, j/2 "C y 3 , 
there can be four more eligible models which are marked by C4 in Table IVl 
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